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We present updated results of the Baksan Underground Scintillator Telescope in search for a signal 
from dark matter particles annihilating in the center of the Sun. Based on the performance of the 
Baksan telescope measuring upward through-going muons with energy larger 1 GeV, we confirm 
the absence of observable excess in arrival directions towards the Sun relatively the expectations 
from neutrinos of atmospheric origin for 24.12 years of live time. We conclude that the 90% 
C.L. upper limits on elastic scattering cross sections of dark matter WIMP (Weakly Interactive 
Massive Particle) off proton are comparable with correspondent limits of other operating neutrino 
telescopes. The best value of the limit is 3 x 10~ 4 picobarn in spin-dependent (SD) interaction of 
WIMP on proton for mass range about 100 200 GeV. 
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1. Introduction 

Neutrino channel in multiwave searches for relic dark matter (DM) particles is accessible 
within data analysis of neutrino telescopes in their regular observations of local sources like the 
Sun where DM could be gravitationally trapped and further accumulated for the time of solar 
system age. Up-to now there are no hints on excess of neutrino events in the direction towards the 
Sun as compared with expected background of atmospheric neutrinos at all neutrino telescopes. 
However, these observations allow to set the upper limits on properties of dark matter particles, in 
particular, the limit on the cross sections of elastic scattering of WIMP on nucleons. Due to the Sun 
chemical composition with approximately 73% of hydrogen, the most sensitive quantity is the spin 
dependent cross sections of dark matter WIMP on proton as it can be seen from the latest results of 
Super Kamiokande [jl|] and IceCube [Q] collaborations. 

The sensitivity of neutrino telescopes strongly depends on their muon energy threshold E'^, 
especially in probing of light masses of WIMP. The Baksan Underground Scintillator Telescope [||] 
(further referred as the Baksan or the BUST) has muon energy threshold around 1 GeV, allowing 
for searches for neutrino signal from annihilation of relatively light WIMPs in the Sun. Minimal 
mass for WIMPs which could be probed at the Baksan is about 10 GeV. This is intrinsic challenge 
in light of recent claimed spectral peculiarities and their interpretations by many authors as an 
evidence on existence of light relic WIMPs both in low background detectors of nuclear recoil 
reactions (see review by J.Cooley at this conference) and in cosmic ray experiments of y, e~e + and 
pp measurements (see review by A.Morselli at this conference). 

Here we present the preliminary results of the Baksan experiment with the statistics twice as 
compared to earlier results [Q], ^ and with a new analysis improved in several ways. Firstly, 
the previous Baksan results P], [j|] were obtained in the framework of a specific model - MSSM 
(Minimal Supersymmetric Standard Model) with neutralino as a dark matter candidate (see review 
of N.Fornengo). Performing the scan over its parameter space but fixing the mass of dark matter 
particle, there have obtained conservative upper limits on muon flux and the rate of dark matter 
annihilation in the Sun. Now rather then studying a specific model we follow present day tendency 
and consider specific dark matter annihilation channels, namely annihilations into bb, and 
VK + IV~. This was also done because the previous tactics makes it difficult to compare the Baksan 
limits with the results of other experiments. Second, we take into account current knowledge about 
neutrino properties in solving the transport of oscillating and interacting neutrinos produced in the 
center of the Sun with energy larger few GeV. 

2. Experiment, triggers and sample of upward through-going muons 

The Baksan Underground Scintillator Telescope is 34-years continuously operating cosmic 
rays detector with 47T-geometry for arrival penetrating charged particles. Main parameters of the 
telescope have been presented in details elsewhere [H], [Qj. Separation of arrival directions between 
up and down hemispheres is made by time-of-flight method with a time resolution 5 ns It was 
shown [[J that 95% of 1 //} values (ratios of speed of light and measured particle track velocity) lie 
in range of 0.7-1.3 ns for single downward going muons. Such interval but with negative sign is 
used to select upward going muons generated by neutrino interactions in down hemisphere. 
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Telescope is located at the altitude 1700 m above see level in the Baksan valley of the North 
Caucause and at the depth 850 g/cm 2 under the mountain Andyrchi where the flux of atmospheric 
downgoing muons is reduced in 5000 times but it is still higher than upgoing muons flux in six 
orders of magnitude. Trajectories of penetrating particles are reconstructed by the positions of hit 
tanks, which put together a system of 3,150 liquid scintillation counters of standard type (70 cm x 
70 cm x 30 cm) in configuration of parallelepiped (17 m x 17 m x 11 m). The counters entirely 
cover all its sides and two horizontal planes inside at the distances 3.6 m and 7.2 m from the bottom. 
Each plane is separated from another one by concrete absorber of 160 g/cm 2 . The configuration 
provides 1.5° of angular accuracy. 

There are two hardware triggers used to select upward going muons. They reduce initial 
rate of downward going muons approximately by factor of 10 3 . Trigger I covers the zenith angle 
range 95° 4- 180° while trigger II selects horizontal muons in the range 80° 4- 100°. The hardware 
trigger efficiency of 99% has been measured with the flux of atmospheric muons. These two 
triggers give about 1,800 events per day for further processing. In the year 2000 the telescope data 
acquisition system was upgraded and allowed to simplify the trigger system down to one general 
trigger with the rate of 17 Hz. All raw information is than undergone further selection using off-line 
reconstruction code. 

Selection criteria of neutrino events in the Baksan neutrino experiment did not changed since 
first data analyses [{||. Few requirements to be satisfied in preliminary selection were determined 
before the telescope was launched: only single trajectory of penetrating particle and a negative 
value of measured j3. Additionally all events with negative values of j8 were scanned by eyes to 
check possible misinterpretation. Also required to have an enter point of each trajectory lower 
than exit point in a range not less than a tank sizes. To reduce background from downward going 
atmospheric muons coming from smaller thickness of mountain but be scattered at large angles, 
all events with azimuths 0° < < 180° in a sample of trigger II are rejected. It was found 1700 
upward through going muons for 24.12 years of live time (l.t.) survived these cuts among detected 
events since December of 1978 till November of 2009. 

Next level of cuts impacts on events mimicked by downward going atmospheric muon inter- 
actions or muon groups: for the reconstructed muon trajectories its path inside the telescope had 
to be larger than thickness 500 g/cm 2 (that corresponds to muon energy threshold about 1 GeV) 
with both enter and exit points; selection only events with —1.3 < 1//3 < —0.7 ns; edge condi- 
tions on enter and exit for events in sample of horizontal events. Totally 1255 events for 24.12 l.t. 
years survived all cuts. From Fig. [I] (left) one can conclude about stable rate of measured upward 
through-going muons during all years of observation. 

3. The Sun survey during three decades 

In search for neutrinos from the dark matter processes in the Sun, we analyse distribution of 
measured events as a function of cos(*F^_sun) where ^-sun is the angle between upward incom- 
ing muons and the Sun position. In Fig. [I] (right) the obtained cosine distribution is presented. The 
mean rate per bin is shown by blue line. We estimate the background expected from atmospheric 
neutrinos directly from real data and shifted (false) Sun positions. Here we follow our previous 
analysis [Q], ^ where it was shown that this method is compatible with Monte Carlo (MC) simu- 
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Figure 1: Left: Measured rate of upward through going muons at the BUST since December 1978. Right: 
Distribution of measured events in angle ^-Sun between incoming neutrino events at the BUST and the 
Sun position. The direction of the Sun corresponds to cos( x I , | 1 _sun) = 1- 



lations of the detector response. Mean energy of simulated neutrinos, which produce muons with 
energy larger 1 GeV to cross the Baksan telescope, is about 50 GeV. From comparison of data 
sample for 21.15 years of l.t. ^ and MC statistics of time larger than real data taking by factor 22 
(i.e., in total of 460 years) it was found that the ratio of observed total number of events to expected 
one without neutrino oscillations is 0.87 ±0.03(stat.) ±0.05(syst.) ±0.\5(theor.). The details has 
been presented elsewhere 

At the location of the Baksan telescope (43,16°N and 42,41°E) the Sun is seen in average 
about half of time per day during a year in both hemispheres; the real bin contents of solar zenith 
(8) distributions reconstructed from data time information (red points with error bars) and one 
(blue histogram) calculated using the Positional Astronomy Library |Q] are shown in Fig. || (left). 
The reconstructed values are obtained from a sample of 1255 selected upward going muons for 
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Figure 2: Left: Zenith distributions of the Sun during the BUST runtime: observed (red points with error 
bars) and calculated (blue, see text). Right: Distributions of cosC^-Sun) for events toward true Sun position 
below horizon (red) in comparison with measured background (green) obtained from the distributions for 
cos( l P^_sun) averaged over six shifted Sun positions. 



4 



Baksan upper limits on SD WIMP-proton cross-sections 



Olga Suvorova 



24.12 years of live time. The calculated distribution is normalized to this number. As it can be seen 
that our multi-years measurements reproduce the Sun full-year passing track with a good accuracy. 
That is the base for our further analysis. 

Potentially the neutrinos from dark matter annihilations in the Sun could be found at night 
time when the Sun is below horizon. For "night" neutrinos we compare cosine distribution in Fig. || 
(right) and integrated angular distribution in Fig. || with respect to measured background which has 




Figure 3: The 90% C.L. upper limits on additional number of events toward the Sun (blue) in cone half- 
angles y. Also shown are measured events (red) and background (green) in the same cones. 

been obtained from averaged distribution of six cases of false Sun positions shifted over ecliptic. 
In Fig. H shown on abscissa value / is a half open angle of the cone toward the Sun. We observe no 
an excess of events coming from the Sun and set the 90% C.L. upper limits on additional events in 
each cone as also shown in Fig. ||. Assuming Poisson statistics for both expected background (Ng) 
and observed events (N b s ) we obtain upper limits C.L. on additional number of signal events (Afc) 
as follows [||] : 

-(Nb+N s ) yKbs (N B +Ns)" 

CL. = 1 X ^ • M 



0-Nb y n °i>* Hsl 
e 2-o 



Further we translate these results into limits on annihilation rates and muon fluxes from the Sun 
direction, considering in details a transport of high energy neutrinos from the Sun to the Baksan 
detector and found dependence on WIMP mass of value y where 90% of N$ events are collected. 



4. Transport of solar neutrinos with energy higher few GeV 

In this Section we briefly describe techniques which we used for numerical simulations of 
muon signal in the BUST and leave the details for future publication [Sj. Resulting neutrino spec- 
trum from WIMP annihilations is a mixture of inclusive spectra from decays of produced copious 
particles. Spectral characteristics of this spectrum depends on the WIMP mass and its annihilation 
channel. To obtain model-independent bounds we consider WIMP of masses from 10 GeV to 1 TeV 
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with three dominant annihilation channels: bb, W + W~ and t~t + . The channel bb represents an 
example of "soft" spectrum while and t + t~ are examples of "hard" spectra. 

We use our own C code for MC simulation of neutrino propagation from the center of the Sun 
to the level of the Baksan detector. At the point of production in the Sun, we take neutrino spectra 
for each annihilation channel obtained in Ref. [|I0|]. When simulating neutrino propagation we take 
into account oscillations in the Sun matter and in the vacuum, CC (charged current) and NC (neutral 
current) scattering of neutrino with interior of the Sun including z regeneration. Our procedure of 



MC simulations is very similar to that of in WimpSIM package [ ITj , H2H. Neutrino oscillations 
in 3 x 3 scheme were implemented according to the algorithm presented in Refs. [13, 14] which 
is very convenient for using with varying density of the Sun. For parameters of neutrino mixing 
matrix we use the best fit of experimental data of Ref. [O]. Below we present the results for normal 
hierarchy only. We use standard solar model presented in Ref. [jl^j. Neutrino-nucleon DIS cross 
sections and distributions are calculated according to formulas presented in Jl7| ] and CTEQ6 JT8| ] 
partem distribution functions. Finally, we calculate flux of muons produced by neutrinos in the rock 



below the Baksan telescope as in Ref. [20]. We check results obtained by our code for neutrino flux 



at the level of the Earth with that of obtained by WimpSIM [ ]1 1 , |12| ] and find an agreement with a 
sufficient accuracy 

5. Results and limits on neutralino-proton cross sections 

If the capture and annihilation of WIMPs in the Sun reach exact equilibrium for time much less 
than age of the solar system 4,5 Giga years, the annihilation rate F A is equal to a half of capture rate 
defined by WIMPs scattering cross sections off solar matter. And in this case, the limit on muon 
flux from dark matter annihilation in the Sun and the limit on annihilation rate can be recalculated 



to limit on elastic cross section of dark matter particle on nucleons pl| , [22] ] . According to that 
annihilation rate can be divided into two pieces 

r A = r A (o S i)+r A {o SD ), (5.1) 

where r^CTs/) and r A (oso) WQ the parts of equilibrium annihilation rate determined by either 
spin-dependent (SD) or spin-independent (SI) interactions. 

To obtain the upper limit on SD elastic cross sections we use recalculation procedure described 



in [22] where we obtained corresponding coefficients in the following expressions 

4T L '' m K) = ^ K) • r U A ppLim (m x ) (5.2) 

(5.3) 

as functions of neutralino mass m x . The (preliminary) results for the limits on SD elastic cross 
section from Baksan's data for the annihilation channels bb and t + t~ are presented in 

Figure |] (left) along with corresponding limits of Super Kamiokande [|l|] and IceCube [Q]. Typically, 
only and bb channels are considered in the limit on dark matter annihilation rate or SD 

elastic cross section. We include also t + t~ channel because of two reasons. Firstly, it allows 
to extend the limits on SD cross section for "hard" type of neutrino spectra to lower masses of 
dark matter. Second, it illustrates an important point: in spite of the fact that the limits on muon 
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Figure 4: The Baksan limits on SD elastic cross section of WIMP on proton in comparison with other 
experimental results (left); ratios of expected muon fluxes both from DM neutrino and antineutrino generated 
in pure three annihilation branches bb quarks, or t + t~ leptons or bosons in cases with and without 

three flavours oscillations (right). 

fluxes are approximately the same for W+W - and t + t~, the limits on SD cross section (and 
annihilation rate) are quite different for these two channels as it can be seen in Fig. | (left). The 
explanation of this fact resides partly in the effect of neutrino oscillations: the enhancement due 
to this effect for t + t~ channel is considerably larger than for the case of VK + W~ channel. This 
is illustrated in Figure |4] (right) where the ratio of muon fluxes with and without oscillation for 
different annihilation channels depending on mass of dark matter is shown. Another important 
effect is a breaking "democracy" in average number of high energy neutrinos per one act of dark 
matter annihilations: for T + T~ pairs it is larger in a factor of 2 4 3 than for or bb channels. 

It can be seen in Fig. |] that neutrino oscillations result in decrease of muon flux from bb 
annihilation channel, which is expected to be dominating at lighter neutralino masses. At the same 
time the muon flux from VK + W~ annihilations is enhanced. These results allow also to interpret 
more correctly the previous Baksan limit [Q] (brown line in Fig. ^, left) calculated without taking 
into account neutrino oscillations. 



6. Summary and conclusions 

We have performed updated analysis of indirect dark matter search with the Baksan Under- 
ground Scintillator Telescope data for 24.12 years of live time. Search for an excess of upward 
going muons in the direction toward the Sun we do not observe any significant deviation from ex- 
pected atmospheric background. We have presented preliminary the 90% C.L. upper limits on spin- 
dependent elastic cross section of dark matter on proton assuming particular annihilation channels 
bb, W + W~ and t + t~. These limits have been derived from the limits on muon fluxes and anni- 
hilation rates for dark matter masses in the interval 104- 1000 GeV. The best value of the limit is 
about 3 x 10~ 4 picobarn for WIMP masses within 1004 200 GeV, that is comparable with those 
presented by Super Kamiokande and IceCube collaborations. 
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